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ABSTRACT 

We present the temperature and polarization angular power spectra of the cosmic 
microwave background (CMB) derived from the first 5 years of WMAP data. The 5- 
year temperature (TT) spectrum is cosmic variance limited up to multipole t = 530, 
and individual ^-modes have S/N > 1 for I < 920. The best fitting six-parameter 
ACDM model has a reduced \ 2 f° r ^ = 33 — 1000 of \ 2 jv = 1.06, with a probability 
to exceed of 9.3%. There is now significantly improved data near the third peak which 
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leads to improved cosmological constraints. The temperature-polarization correlation 
(TE) is seen with high significance. After accounting for foreground emission, the \ow-£ 
reionization feature in the EE power spectrum is preferred by Ax 2 = 19.6 for optical 
depth r = 0.089 by the EE data alone, and is now largely cosmic variance limited 
for £ = 2 — 6. There is no evidence for cosmic signal in the BB, TB, or EB spectra 
after accounting for foreground emission. We find that, when averaged over t = 2 — 6, 
£{£ + l)Cf B /{2ir) < 0.15 /iK 2 (95% CL). 

Subject headings: cosmic microwave background, cosmological parameters, cosmology: 
observations, early universe, large-scale structure of universe, space vehicles: instru- 
ments 



Introduction 



The WMAP satellite (IBennett et al.l 120031 ) has measure d the temperature an d polarization 
of the microwave sky at five frequencies from 23 to 94 GHz. iHinshaw et al.l (J2008J) presents our 
new, more sensitive tempe rature and polariz ation maps. After removing a model of the foreground 
emission from these maps (jGold et al.ll2008l ). we obtain our best estimates of the temperature and 
polarization angular power spectra of the cosmic microwave background (CMB). 

This paper presents our statistical analysis of these CMB temperature and polarization maps. 
Our basic analysis approach is similar to the approach described in our first year W MAP tem- 
perature analysis ( Hinshaw et al. 20031') and polariza tion analysis ( Kogut et al. 20031') and in th e 
three year WMAP temperature (jHinshaw et al.ll2007l ) and polarization analysis (jPage et al.ll2007l ). 
While most of the WMAP analysis pipeline has been unchanged from our 3-year analysis, there 
have been a number of improvements that have redu c ed th e systematic errors and increased the 
precision of the derived power spectra. Hinshaw et al] ( 2008 ) describe the WMAP data processing 
with an emphasis on these changes. iHill et al.1 (|2008l ) present our more complete analysis of the 
WMAP beams based on 5 years of Jupiter data and physical optics fits to both the A- and B-side 
mirror distortions. The increase in main beam solid angle leads to a revision in the beam func- 
tion that im pacts our compute d power spectrum by raising the overall amplitude for £ > 200 by 
roughly 2%. iGold et al.1 (120081 ) introduce a new set of masks that are designed to remove regions 
of free-free emission that were a minor (but d etectable) contamin ant in analyses using the previous 
Kp2 mask u sed in the 1- and 3-yea r analysis. IWright et al.l (|2008l ) updates the point source catalog 
presented in IHinshaw et al.l (|2007l ) finding 67 additional sources. Section [2] of this paper describes 
these changes and their implications for the measured power spectra. 

Section [3] presents the temperature angular power spectrum (TT). WMAP has made a cosmic 
variance limited measurement of the angular power spectrum to £ = 530 and we now report results 
into the "third peak" region. The WMA P results, combined with recent ground-based measure 



ments of the TT angular power spectrum ([Readhead et al.ll2004l ; iJones et al.ll2006l ; iReichardt et al 
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2008), result in accurate measurements well into the "fifth peak" region. For WMAP, point sources 
are the largest astrophysical contaminant to the temperature power spectrum. We present estimates 
for the point source contamination based on multi-frequency data, source counts and estimates from 
the bispectrum. 



The polarization observation s are decomposed into E and B mode components (jKamionkowski et al 



19971 : ISeliak fc Zaldarriagalll997l ). Primordial scalar fluctuations generate only E modes, while ten- 
sor fluctuations generate both E and B modes. With T, E and B maps, we compute the angular 
auto-power spectra of the three fields, XT, EE and BB, and the angular cross-power spectra 
of these three fields, TE, TB and EB. If the CMB fluctuations are Gaussian random fields, 
then these six angular power spectra encode all of the statistical information in the CMB. Un- 
less there is a preferred sense of rotation in the universe, symmetry implies that the TB and EB 
power spectrum are zero. In Section [4] we present both the TE and TB temperature-polarization 
cross power spectra. The WMAP measurements of the TE spectrum now clearly see multiple 
peaks. The large angle TE a n ti-correlation is a dist inctive signature of superhorizon fluctuations 
(jSpergel fc Zaldarriagalll997al ). iKomatsu et al.l (120081 ) discuss how the TB measurements constrain 
parity- violating interactions. Section [5] presents both the EE and BB polarization power spectra. 
The E E power spectrum now shows a clear ~ 5a signature of cosmic reionization. iDunklev et al 



([20081) s how that the amplitud e of the signal implies that the cosmic reionization was an extended 



process. 

the angular power spectrum measurements. 



Dunklev et al.l (j2008l ) and IKomatsu et al.l (j2008l ) discuss the cosmological implications of 



2. Changes in the 5-year Analysis 

The methodology used for the 5-year power spectra analysis is similar to as that used for the 
3-year analysis. In this section we list the significant changes and their impact on the results: 



• iHinshaw et al.l (j2008l ) describe the changes in the map processing and the resultant reduction 
in the absolute calibration uncertainty from 0.5% to 0.2%. 



The temperature mask used to compute the power spectrum has been updated 



slight ly more sky near the galactic plane, and more high-latitude point sources (jGold et al 



removing 



20081 ). The galactic mask used in the 3-year and 1-year releases (Kp2) was constructed by 
selecting all pixels whose K-band emission exceeded a certain threshold. This procedure 
worked well in identifying areas contaminated by synchrotron emission; however, it missed 
a few small regions contaminated by free-free, particularly around p Oph, the Gum nebula, 
and the Orion/Eridanus Bubble. For the 5-year analysis we have constructed a new galactic 
mask to remove these contaminated areas. 



Wright et al.l (|2008l ) updated the WMAP point source catalog, finding 390 sources in the 
5-year data, 67 more sources than in the 3-year catalog (jHinshaw et al.l 120071 ). Of these 67 
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new sources, 32 were previously unmasked, and therefore added to the 3-year source mask to 
create the 5-year source mask^j 

All told, the new 5-year temperature power spectrum mask (KQ85) retains 81.7% of the sky, 
while the 3-year mask (Kp2) retained 84.6% of the sky. 



• The 5 -year polarization mask is the same as the 3-year P06 mask described in iPage et al 
([20071 ) . except that an additio nal 0.27% of the sky has been removed due to combining P06 
with the new processing mask (IHinshaw et al.ll2008l ). 



In addition to masking, the maps are further cleaned of galactic foreground emission using 
ext ernal templ a tes. T he cleaning procedure is very similar to that of the 3-year analysis; 



sec 



Gold et al.l (|2008l ) for details. For the temperature map, three templates are used: a 
synchrotr on template (the WMAP K — Ka difference map ), an Ha template as a proxy for 



free- free (iFinkbeinerl 120031 ) . and a thermal dust template (IFinkbeiner et al.lll999l ). For the 
polarization maps, two templates are used (since free- free is unpolarized): the polarized ri- 
band map and a polarized dust template constructed from the unpolarized dust template, a 
simple model of the galactic magnetic field, and polarization directions deduced from starlight. 

A great de al of work has go ne into improving the determination of the beam maps and window 
functions (|Hill et al.l l2008). The main beam solid angles are larger than the 3-year estimates 
by as 1-2% in V- and W-band. Increased solid angle (i.e., greater map smoothing) reduces 
the value of the transfer function bi, raising the deconvolved CMB po wer spectra . The ratio 
of the 3-year to 5-year transfer functions can be seen in Figure 13 of iHill et al.1 (|2008l ): the 
net effect is to raise the TT power spectrum by ~ 2% for I > 200, which is within the 3-year 
beam la confidence limits. The beam transfer function uncertainty is smaller than the 3-year 
uncertainty by a factor of ~ 2. The window function uncertainty is now ~ 0.6% in AC//C/ 
for 200 < t < 1000. 



3. Temperature Spectrum 



The 5-year I < 32 spectrum is described in iDunkley et al.1 (120081 ) . At low-£ the likelihood 
function is no longer well approximated by a Gaussian so that we explicitly sample the likelihood 
function to evaluate the statistical distribution of each multipole. 

We construct the 5-year TT spe c trum for I > 32 in the same fashion as the 3-year spectrum; 
we refer the reader to IHinshaw et al.l (|2007l ) for details, and only briefly summarize the process as 
follows: 



Six of the sources in the 5-year catalog were not added to the mask; they were found in a late update to the 
catalog after the mask had been finalized. 
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• We start with the single year V1,V2,W1-W4 resolution-10 mapsp masked by the KQ85 mask, 
and further cleaned via foreground template subtraction. 

• The pseudo-C; cross power spectra are computed for each pair of maps. Two weightings are 
used: flat weighting and inverse noise variance (iV b s ) weighting. 

• The year/DA cross power spectra are combined by band, forming the Vx V, VxW, and Wx W 
spectra. The auto power spectra are not included in the combination, eliminating the need 
to subtract a noise bias. 

• A model of the unresolved point source contamination with amplitude A ps = 0.011±0.001 /xK 2 sr 
is subtracted from the band-combined spectra. See Section [37X1 for more details. 

• The VxV, VxW, & WxW spectra are optimally combined £-by-£ to create the final CMB 
spectrum. 



As in the 3-year analysis, the diagonal elements of the C\ covariance matrix are calculated as 

where C\ is the cosmic variance term and Ni the noise term. The value of / s ky(0> the effective sky 
fraction, is calibrated from simulations! 

f 0.826 - 0.091(^/500) 2 , £ < 500; 
/sky( ) - I Q 7?7 _ .i27(500/£), £ > 500. ^ ' 

The 5-year TT spectrum is shown in Figure [TJ With the greater S/N of the 5-year data the 
third acoustic peak is beginning to appear in the spectrum. The spectrum is cosmic variance limited 
up to £ = 530, and individual £- modes have S/N > 1 for £ < 920. In a fit to the best cosmological 
ACDM model, the reduced x 2 f° r ^ = 33 — 1000 is x 2 / u = 1-06, with a probability to exceed of 
9.3%. 

Figure fl2l compares the unbinned 5-year TT spectrum with the 3-year result. Aside from the 
small upward shift of the 5-year spectrum relative to that of the 3-year, due to the new beam transfer 
function, they are identical at low-£. Figure [TBI shows the unbinned TT spectrum broken down into 
its frequency components (VxV, VxW, WxW), demonstrating that the signal is independent of 
frequency. 



2 12,582,912 pixels (7V sidc = 1024). 

3 The Markov chains in Dunklev et al. 1 2008h and Komatsu et al. (2008?) were run with a version of the WMAP 
likelihood code with older and slightly larger values for / s k y . The change in / s it y increased the TT errors by on average 
2%. Rerunning the ACDM chain with the new / s k y leads to parameter shifts of at most O.lo". 
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How much has the determination of the 3rd acoustic peak improved with the 5-year data? 
Over the range I = 680 — 900, which approximately spans the rise and fall of the 3rd peak (from 
the bottom of the 2nd trough to the point on the opposite side of the peak), the fiducial spectrum 
is preferred over a flat mean spectrum by A% 2 = 7.6. For the 3-year data it was A% 2 = 3.6. With 
a few more years of data, WMAP should detect the curvature of 3rd peak to greater than 3a. 

In Figure [2] we com pare the WMAP 5-y e ar TT power spectrum along with rece nt results 
from other experiments (jReadhead et al.l [200J; iJones et al.l 120061 ; iReichardt et al.l 120081 ) , showing 
great consistency between the various measurements. Several on-going and future ground-based 
CMB experiments plan on calibrating themselves off their overlap with WMAP at the highest-£'s; 
improving WMAP's determination of the 3rd peak will have the added benefit of improving their 
calibrations. 



3.1. Unresolved Point Source Correction 

A population of point sources, Poisson-distributed over the sky, contributes an additional 
source of white noise to the measured TT power spectrum, Cj T — > Cj T + C ps . Given a known 
source distribution N(> S), the number of sources per steradian with flux greater than S, the 
point-source induced signal is 

fSc rfAT 

C ps = g{vf \ dS ^S 2 [/zK 2 sr] (3) 



Jo dS 

where S is the source flux, S c is the flux cutoff (above which sources are masked and removed from 
the map), and g(v) = (c 2 /2ku 2 )r(v) converts flux density to thermodynamic temperature, with 

(e x - l) 2 

r (^) = x 2 e x — ' x = W^cmb (4) 
converting antenna to thermodynamic temperature. 

At the frequencies and flux densities relevant for WMAP, source counts are dominated by flat- 
spectrum rad io sources, wh i ch ha ve flux spectra that are nearly constant with frequency (S ~ u Q 



with a ~ 0). IWright et al.l (|2008l ) finds the average spectral index of sources bright enough to be 
detected in the WMAP 5-year data to be (a) = —0.09, with an intrinsic dispersion of a a = 0.176. 
Since a source with flux S ~ v a has a thermodynamic temperature T ~ v a ~ 2 r(v), we model the 
frequency dependence of C ps as 

C ps {u u uj) = A ps r(v t )r(vi) j (5) 

where Vij are the frequencies of the two maps used to calculate the TT spectrum, A ps is an unknown 
amplitude, and vq = 40.7 GHz is the Q-band central frequency. 



In this section, we estimate the value of A ps needed to correct the TT power spectrum, finding 

ps 



Aps = 0.011 ± 0.001 ^K 2 sr, and discuss incorporating its uncertainty into the likelihood function 
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3.1.1. Estimating the correction 



For a fixed beam size, flat-spectrum radio sources are much fainter in the W-band temperature 
maps than in Q- or V-band, allowing us to use the frequency dependence of the TT spectrum at 
high-£ to constrain the value of A ps . As in previous releases, the estimator we use is 



A 



ps 



hi 



n 



-Ua/3 



(6) 

(7) 



where greek letters represent a pair of frequencies (e.g., VW), CP is the measured TT cross-power 



spectrum, is the (C"C ; P ) covariance matrix including cosmic variance and detector noise, and 
sf = 1(1 + l)C ps (a)/2ir. The inverse estimator variance ([<L4 ps ]~ 2 ) is given by the denominator of 
([6]). While S"' 3 does not include the off-diagonal coupling due to the mask, the diagonal elements 
are renormalized to account for the loss of sky coverage. 

Measured values for A ps are listed in Tabled] for various frequency combinations (QVW & VW) 
and galactic masks (KQ85, KQ80, & KQ75). The QVW estimates are insensitive to the galactic 
mask; the VW estimate increases somewhat as more of the sky is masked. Both the QVW and VW 
estimates prefer the same value (~ 0.011 /iK 2 sr) of A ps when the KQ75 mask is used. While we 
restrict the data to i = 300 — 800, the QVW estimate is only a weak function of the chosen ^-range; 
Figure [3] shows A ps estimated in bins of width A£ = 100. We adopt A ps = 0.011 ± 0.001 /iK 2 sr as 
our correction to the final combined TT spectrum. The consistency between £-b'ms and between 
QVW and VW seen in Figure[3]is an important null test for the angular power spectrum. A ps (VW) 
is proportional to the power in the (V-W) map in a given Grange. Figure 0] shows no evidence for 
any detectable residual signal in the VW maps after point source subtraction. 

Because radio sources can only have positive flux they intr oduce a positiv e skew ness to the 
maps, which can be detected in searches for non-Gaussianity. iKomatsu et al.l (|2008l ) estimated 
the bispectrum induced by sources, finding V s = (4.3 ± 1.3) x 10" 5 /xK 3 sr 2 at Q-band. Is this 
consistent with the value of A ps measured from the power spectrum? Given a theoretical model 
for the source number counts iV(> S), one can predict the meas ured values of C ps and Several 
models exist in the li terature; we tested our results against two, Toffolatti et al. (1998, Tof98)0 and 
de Zotti et all iboOfJ . deZ05). C ps is calculated via and V s from 







(8) 



Bennett et al. (2003) w e found that the Tof98 model needed to be rescaled by a factor of 0.66 to match the 



WMAP 1-year number counts; Wright et al. ( 20081 ) refined the rescaling factor to 0.64 to match the WMAP 5-year 
source counts. 
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where g{y) and S c are denned in ([3]). The comparison is complicated by the fact that S c is 
unknown. We mask out not only the sources detected in WMAP data, but also undetected sources 
from external catalogues that are likely to contribute contaminating flux. However, a single value 
of S c predicts both C ps and W s , so we can in principle tune S c to match one, and see if it agrees 
with the other. In Table [2] we compare our measured values of C ps and with the rescaled Tof98 
and deZ05 predictions for several values of S c . There is some tension between the measured values 
and the model predictions. Given our measured value for b ps the models would prefer a smaller 
value for A ps , in the range 0.008-0.010 /iK 2 sr. For the Tof98 model, the S c 0.52 predictions are 
within la of both C ps and V s . However, the deZ05 model appears to be discrepant, and a single 
value for S c cannot match both C ps and bP s . 

Other groups have independently estimated the unresolved source contamination, and their 
results are in general agreement with ours. When the 3-y e ar da ta was initially released the correc- 



tion was A ps = 0.017 ± 0.002 /xK 2 sr. iHuffenberger et al.l ( 2006 ) reanalyzed the data and claimed 



A ps = 0.011 ± 0.001, noticing that A ps was sensitive to the choice of galaxy mask; using the KpO 
mask instead of Kp2 reduced the value of A ps . Revisiting our original estimate for the 3-year anal- 
ysis, we reduced the corr ection to 0.014 ± 0.003 for the published papers. In a subsequent paper, 



Huffenberger et al.l (|2007l ). the same group corrected their original estimate after finding a small 



error, finding 0.013 ± 0.001, consistent with our published result. 



4. Temperature-Polarization Spectra 



The standard model of adiabatic primordial density fluctuations predicts a correlation between 
the temperature and polarization fluctuations. The temperature traces primarily the density, and 
E-mode polarization the velocity, of the photon-baryo n pl asma at rec o mbin ation. The correlation 
was seen in earlier WMAP data by lKogut et al.l (|2003l ) and lPage et al.l (120071 ). The anti-correlation 
near t = 30 provides evidence that fluctuations exist on super horizon scales, as it is o bserve d on 
an angular scale larger than the acoustic horizon at decoupling ISpergel Zaldarriagal (|1997bl ). 



No significant c hanges have been made in the five-year TE analysis. We continue to use the 
method described in iPage et al.1 (1200711 to compute the TE power spectrum. The inputs are the 
KaQV polarization maps (|Gold et al.ll2008l ). and the VW temperature maps. For high multipoles 
£ > 23, the likelih ood can be approx imated as a Gaussian, and we continue to use the ansatz given 
in Appendix C of lPage et al.1 (120071 ) to compute the covariance matrix. At low multipoles, i < 23, 
th e likelihood of the polarization data is evaluated directly from the maps, following Appendix D 
in 



Page et al.l (120071 1. 



Figure [5] shows the TE spectrum. At low-£ the spectrum and error bars are approximated using 
the Gaussian form, although these are not used for cosmological analysis. With five years of data 
the anti-correlation at t = 140 is clearly seen in the data, and the correlation at i = 300 is measured 
with higher accuracy. The second anti-correlation at I ~ 450 is now better characterized, and is 
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consistent with predictions of the ACDM model. The structure tests the consistency of the simple 
model, which fits both the TT and TE spectra with only six parameters. The best-fit ACDM 
model has x 2 = 415 for the TE component, with 421 degrees of freedom, giving x 2 l v = 0.99. 
The consistency confirms that the fluctuations are predominantly adiabatic, and constrains the 
amplitude of isocurvature modes. 

The signal at the lowest multipoles, evaluated using the exact likelihood, is used to provide 
additional constraints on the reionization history. Although small, the measur ement is co n sisten t 
with the EE s i gnal, and consistent with the three- year WMAP observations iPage et al.l (|2007l ). 
Dunklev et al.l (|2008l ) discuss constraints on reionization. 



No correlation is expected between the temperature and the B-mode polarization. The TB 
spectrum is therefore primarily used as a null test, and is shown in Figure El It is consistent with 
no signal, as expected ; over £ = 24 — 450 the reduced null x 2 is 0.97. This measurement is used in 
Komatsu et al.l (|2008l ) to place constraints on the presence of any parity violating terms coupled to 
photons, that could produce a TB correlation. We now include the TB spectrum at high £ as an 
optional module for the likelihood code. 



5. Polarization Spectra 



Due to its thermal stability (jjarosik et al.1 120071 ) and well-characterized gain, WMAP can 
measure polarization signals even though the scan pattern was not optimized for doing so. The 
polarization signal is manifested in the time ordered data (TOD) differently from the temperature 
signal. As a result, some of the low-£ polarization mul tipoles are well sampled and other mu ltipoles 
are poorly sampled and have large statistical errors (jHinshaw et al.l 120071 ; IPage et al.l 120071 ) . This 
is a rather different situation than from that of the temperature spectrum, and the data must be 
analyzed with some care. 

When we analyze the I = 2 temperature power spectrum, we use the likelihood function rather 
than Gaussian errors, as the Gaussian approximation starts to break down with only ~ 4 effective 
modes measured in the map (the reduction is due to / s k y ~ 0.7). For polarization, this effect is even 
more dramatic, as our scan pattern significantly lowers the effective number of multipole s measured 
parti cularly for EE £ = 2,5,7 and 9 and BB 1 = 3 (the peaks seen in Figure 16 in IPage et al 



20071 ). Figure [8] demonstrates the importance of using the full likelihood description. The Figure 
shows both the pseudo-C/ estimates of the I = 2 — 7 BB multipoles and the conditional likelihoods 
computed using the WMAP likelihood code by varying the multipole in question, keeping the rest 
of the spectrum fixed to the fiducial best-fit ACDM model. From the plots it is clear that the best 
estimates of the mean and the uncertainty are not attained with the pseudo-C; estimates. 

We next consider the low-£ EE and BB power spectra in more detail. The low-£ EE power 
spectrum is shown in Figure [9j The uncertainties are obtained from the conditional likelihood and 
include cosmic variance; thus one cannot double the error flags to get the 95% confidence limits. 
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If we zero out the £ < 10 portion of the fiducial EE & TE spectra the x 2 increases by 22.3, of 
which 2.7 is due to TE. Thus the reionization feature in the EE power spectrum is preferred by 
A X 2 = 19.6. The I = 2, 3, 4, & 6 multipoles are cosmic variance limited, and the S/N ratio for the 
combined £ = 2 — 7 bandpower is 11. 

Considerable effort has gone into understanding the W-band £ = 7 EE signal. Because of the 
apparent anomalously high £ = 7 EE value computed by the pseudo-C^ algorithm, we have avoided 
using the W-band maps in cosmological analysi s and use them only as an additional check on 
various models. Figure 8 of iHinshaw et alJ ((2008) shows that the £ = 7 value, while high, appears 
to be consistent with being in the tail of a properly computed likelihood distribution. The W-band 
£ = 7 problem may be a signature of poor statistics rather than a systematic. However, more data 
are needed to understand this potential anomaly. The £ = 3 BB signal gives perhaps the clearest 
example of the importance of using the full likelihood code. While the pseudo-Q estimate implies 
a significant detection of power, the full likelihood code shows this to not be the case. The physical 
cause of the large uncertainty is that with our scan strategy an £ = 3 BB signal resembles an offse t 



in the data and thus is not well separated from the baseline (jPage et al 



2007 



Hinshaw et al.ll2008h . 



We see no evidence for a B-mode signal at low £, limiting the possible level to ^(^+1)C^^_ 6 / (2ir) < 
0.15 /j,K 2 (95% CL), including cosmic variance. With r = 0.1 and r = 0.2, a typical estimate for 
currently favored models of inflation, 1(1 + l)CfJ 2 _ 6 /(2n) ps 0.008 fiK 2 . Since a signal of 0.15 fxK 
corresponds roughly to r ~ 20, one can see that WMAP's limit is n ot based on the BB da ta, but 



on the tensor contribution to the TT and EE spectra as discussed in lKomatsu et alJ (|2008l ). 



For EE at £ > 10, there are hints of signal in the data consistent with the standard ACDM 
model. However, the significance is not great enough to contribute to knowledge of the cosmological 
parameters. The 5-year high-i 1 EE spectrum is shown in Figure 1 101 along with recent r esults 



from ground-based experiments (jLeitch et al 



2005 



Montroy et al. 



Sievers et al.ll2007h . For 



£ = 50-800, x 2 = 859.1 assuming Cf E = 0, and drops by 8.4, or almost 3cr, assuming the standard 
ACDM model. For the 3-year data the equivalent change in x 2 was 6.2. 

The high-£ BB spectrum is consistent with no signal, having a reduced x 2 °f 1-02 over £ = 
50 — 800 for the QV data. The lack of any signal in the low and high £ BB data is a necessary 
check of the foreground subtraction. As seen in Page et al. (2007), foreground emission produces 
E-modes and B-modes at similar levels; thus the absence of a B-mo de signal suggests tha t the level 
of contamination in the E-mode signal is low. This is quantified in lDunkley et alJ (|2008l ). 



6. Summary and Conclusions 

We have presented the temperature and polarization angular power spectra of the cosmic mi- 
crowave background (CMB) derived from the first 5 years of WMAP data. With greater integration 
time our determination of the third acoustic peak in the TT spectrum has improved. The low-£ 
reionization feature in the EE spectrum is now detected at nearly 5a. The TB, EB, &: BB spectra 
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show no evidence for cosmological signal. The spectra are in excellent agreement with the best 
fit ACDM model. Our knowledge of the power spectrum is improving both due to more detailed 
analyses, better modeling and understanding of the foreground emission, and more integration 
time. 

All of the 5-year WMAP data products are being made available through the Legacy Archive 
for Microwave Background Data Analysis (LAMBDA^), NASA's CMB Thematic Data Center. The 
temperature and polarization angular power spectra presented here are available, as is the WMAP 
likelihood code which incorporates our estimates of the Fisher matrix, point sources and beam 
uncertainties. 



The WMAP mission is made possible by the support of the Science Mission Directorate Office 
at NASA Headquarters. This research was additionally supported by NASA grants NNG05GE76G, 
NNX07AL75G SOI, LTSA03-000-0090, ATPNNG04GK55G, and ADP03-0000-092. EK acknowl- 
edges support from an Alfred P. Sloan Research Fellowship. This research has made use of NASA's 
Astrophysics Data System Bibliographic Servi ces. We acknowledge use of the CAMB, CMBFAST, 
CosmoMC, and HEALPix (jGorski et al.ll2005l ) software packages. 



A. Likelihood treatment of source/beam uncertainties 

In this section, we test the treatment of the unresolved source correction and beam uncertainties 
in the WMAP likelihood code, and show that it produces the correct results for cosmological 
parameters. 

We adopt the same likelihood treatment of the un resolved point source correction uncertainty 



for the 5-year likelihood code as used in the 3-year code (jHinshaw et al.ll2007l . Appendix A), updated 
for the 5-year value of A ps . Briefly, a correction to the logarithmic likelihood, C = — 21nL = Co+C\ 
where Co is the standard likelihood and C\ the combined source & beam correction, is calculated 
assuming the Ci are normally distributed, a reasonable assumption at high I. 



Huffenberger et al.1 (120071 . Huf08) disagreed with the source & beam likelihood module used 
in the 3-year analysis, pointing out that the uncertainty in n s (the index of primordial scalar 
perturbations) was unchanged even if the uncertainty in A ps was increased by a factor of 100 
(Fig. 2 in their paper). They proposed an alternative approach, integrating the beam/point source 
covariance matrix into the cosmic variance/noise/mask covariance matrix and inverting the result 
in order to compute C directly, instead of calculating L\ as a separate correction. Using this form of 
the likelihood, as 5A ps was increased, the uncertainty in n s increases (albeit modestly; 5n s increased 
by 38% when 5A ps -» 100 x 5A ps ). 



http : //lambda. gsfc .nasa.gov/ 
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However, while we agree that it is striking that the error in n s is seemingly unaffected by the 
uncertainty in A ps , we have some concerns regarding the Huf08 approach. To quote Huf08, "the 
errors on the source measurement do not make much difference, as long as [<L4 ps ] < 0.003 [ / uK 2 sr]", 
and their Fig. 2 implies the same holds true when SA ps = 0.003. This value is significant, because 
it is the uncertainty adopted for the 3-year WMAP analysis. When Huf08 adopted the same 
uncertainty, they found the same absolute uncertainty in n s as the WMAP team, but their central 
value was shifted higher by 0.005. This shift persisted as 5A ps — ► 0, and thus was seemingly not 
due to the point source uncertainty. The conclusion we draw is that they found the same value of 
5n s as the WMAP 3-year analysis, but their value of n s was biased high because of the way they 
treat the beam uncertainties. We believe the Huf08 value of n s would be in agreement with that 
found in WMAP3, but that it is biased high due to their treatment of beam uncertainties. 

Huf08 quoted the value of C\ computed with their alternative likelihood module for a particular 
CMB spectrum distributed with the WMAP 3-year likelihood code test program, finding C\ = 
—2.64, whereas the WMAP value is C\ = —1.22. As a check, we numerically marginalize the Cq 
portion of the likelihood over beam and point source errors, to see if we can reproduce their value. 
The desired integral is 

exp(-A/2) = J dxdye-^+y^LoidlQ&y)) (Al) 

where 

Q(x, y) ^ (cf T + xaf tsrc ) (l + £ yi a\^{i) j , (A2) 

is the theoretical model (C, ) perturbed by point source and beam errors. With 10 dimensions 
to integrate over (nine beam modes and one point source mode), normal grid-based quadrature is 
impractical, so we turn to Monte Carlo integration instead: 

-, N M c 

e X p(-£^ C /2) « — Y e m^(<i|C i (xW,^)))-lnL( d |C ( ) (A3) 



Nmc . . 
t=i 



where and yp are independent unit-variance normal deviates. With Nmc = 10 4 points, we 
find Cf c = -1.29 ± 0.04, consistent with the WMAP result of -1.22, but not the Huf08 result of 
-2.64. 

As a further test of whether the our cosmological parameter estimates fully capture the point 
source uncertainty, we have run a Markov chain with a modified form of the point source likelihood 
module, dubbed SRCMARG. The point source correction is calculated via a simple numerical 
integration, 



exp(-£f sr 72) = [ da-^e-^^LoidlCt+aa? 3 ™) 

J v27T 



(A4) 



A N 

~ w ie -^ 2 / 2 L (d\C l + iAaf tsrc ) (A5) 



2 ^ N 
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with N = 25, A = 0.2, and u>j = 1 for except at the endpoints where w\m = 1/2 (the trapezoidal 
rule). The resulting one-dimensional marginalized distribution for erg, shown in the left panel of 
Figure [TT1 is indistinguishable from our standard result. We have also run a SRCMARG chain with 
the error increased by a factor of 5 (i.e., 5A ps = 0.005). In this case the uncertainty in erg increases 
by 15%. 

Likewise, we have run similar tests of the beam uncertainty, dubbed BEAMMARG. The ap- 
proach is the same as SRCMARG, but with "C, + aaf tsrc " in (JAU) re placed by "G(l + aof eam )", 



where af eam is the noisiest beam eigenmode, shown in Figure 12 of iHill et al.l ( 20081 ). The ID 
marginalized distributions for n s are shown in the right panel of Figure [TTJ As with SRCMARG, 
the BEAMMARG result is indistinguishable from our standard result. Inflating the beam error by 
a factor of 20 results in a 14% increase in Sn s , along with a slight shift in n s away from unity. 
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Table 1. Unresolved Point Source Contamination 



Bands 


Mask 


A ps (a = 0) [10- 3 /uK 2 sr] 


A ps (a = -0.09) [10- 3 /uK 2 sr] 


QVW 


KQ85 


11.3 ±0.9 


11.2 ±0.9 




KQ80 


11.3 ±0.9 


11.2 ±0.9 




KQ75 


10.7 ± 1.0 


10.6 ± 1.0 


vw 


KQ85 


6.9 ±3.4 


7.2 ±3.5 




KQ80 


9.1 ±3.6 


9.5 ±3.8 




KQ75 


10.5 ±3.9 


11.1 ±4.1 



Note. — All results are for I = 300 - 800. 



Table 2. Unresolved Point Source Contamination 





S c [Jy] 


C ps [10- 3 /iK 2 sr] 


W s [10~VK 3 sr 2 ] 


WMAP5 (KQ75) 




11.7 =b l.l a 


4.3±1.3 b 


Toffolatti et al. (1998^1x0.64 


0.6 


12.1 


6.8 




0.5 


10.4 


4.9 


de Zotti et al. (2005) 


0.7 


11.7 


8.4 




0.5 


8.3 


4.3 



a By equation ([5]), C ps (Q) = A ps r{Q) 2 = 1.089 x A ps , where A ps is the 
QVW/KQ75 result from Tabled) 



b From iKomatsu et al.1 ( 20081 ). using the Q-band map and KQ75 mask. 



Note. — All numbers are evaluated at 40.7 GHz (Q band). 
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Table 3. Beam/source likelihood treatment effect on parameters 



Likelihood treatment 


n s 






08 




standard 


0.964 


± 


0.014 


0.796 


± 


0.036 


SRCMARG 


0.964 


± 


0.014 


0.798 


± 


0.036 


SRCMARG x5 


0.965 


± 


0.015 


0.803 


± 


0.042 


BEAMMARG 


0.964 


± 


0.015 


0.799 


± 


0.036 


BEAMMARG x20 


0.958 


± 


0.016 


0.796 


± 


0.034 



Note. — One-dimensional marginalized values for n s 
and as for various treatments of the unresolved point 
source and beam uncertainty in the WMAP likelihood 
code. See Section [A] for descriptions of SRCMARG and 
BEAMMARG. Here "x5" and "x20" indicate the error 
has been increased by a factor of 5 and 20, respectively. 
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Fig. 1. — The WMAP 5-year temperature (TT) power spectrum. T he red curve is the best-fit 
theory spectrum from the ACDM/ WMAP chain (jPunklev et alJbood . Table 2) based on WMAP 
alone, with parameters {VL b h 2 , n m h 2 , A|,, n s , r, H ) = (0.0227,0.131,2.41,0.961,0.089,72.4). The 
uncertainties include both cosmic variance, which dominates below £ = 540, and instrumental 
noise which dominates at higher multipoles. The uncertainties increase at large I due to WMAP's 
finite resolution. The improved resolution of the third peak near I = 800 in combination with the 
simultaneous measurement of the rest of the spectrum leads to the improved results reported in 
this release. 
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Fig. 2. — The WMAP 5-year TT power sp ectrum along wit h recent results fr om the ACBAR 
(Rei chardt et alj|2008l . purple), Boomerang (| Jones et al.ll2006l . green), and CBI (jReadhead et al 



2004j , red) experiments. The other experiments calibrate with WMAP or WMAP's measurement 
of Jupiter (CBI). The red curve is the best-fit ACDM model to the WMAP data, which agrees well 
with all data sets when extrapolated to higher-^. 
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Fig. 3. — The unresolved point source contamination A ps , measured in bins of Al = 100 evaluated 
at 40.7 GHz (Q-band). For a source population whose fluxes are independent of frequency A ps 
scales roughly as ~ v~ 2 in the WMAP data. The red data points are from the analysis of V and 
W bands alone and the blue points are from the analysis of Q, V, and W bands. The horizontal 
dashed green lines, at 0.010 and 0.012, show the la bounds for our adopted value of A ps . Note 
that the QVW amplitude is independent of £. 
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Fig. 4. — The TT V — W null spectrum. After correcting for unresolved point source emission, 
the individual power spectra are subtracted in power spectrum space. The result is consistent with 
zero and thus there is no evidence of point source contamination. In these units, point source 
contamination would be evident as a horizontal offset from zero. At £ = 500, the TT power 
spectrum is Cj T ~ 0.06; thus the contamination is limited to roughly 3% in power. 
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Fig. 5. — The WMAP 5-year TE power s pectrum. The green curve is the best-fit theory spectrum 
from the ACDM/WMAP Markov chain (jPunkley et alJhood ). For the TE component of the fit, 
X 2 = 415, and there are 427 multipoles and 6 parameters; thus the number of degrees of freedom 
is v = 421, leading to x 2 l v = 0.99. The particle horizon size at decoupling corresponds to / ~ 100. 
The clear anticorrelation between the primordial plasma density (corresponding approximately to 
T) and velocity (corresponding approximately to E) in causally disconnected regions of the sky 
indicates that the primordial perturbations must have been on a superhorizon scale. Note that the 
vertical axis is (£ + l)Q/(2vr), and not 1(1 + l)Q/(2vr). 



- 22 - 



2.0 


1 1 1 


1.5 




1.0 




0.5 


< > 







0.5 




1.0 


1 1 1 



"i — i — i — r 



i 1 r 



j i 



J I I I I I I L 



10 



100 500 

Multipole moment / 



1000 



Fig. 6. — The WMAP 5-year TB power spectrum, showing no evidence of cosmological signal. The 
null reduced \ 2 f° r ^ = 24 - 450 is 0.97. Note that the vertical axis is (£ + l)Q/(27r), and not 
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Fig. 7. — Conditional likelihoods for the t = 2 — 7 EE multipole moments (black curves), com- 
puted using the WMAP likelihood code by varying the multipole in question, with all other 
multipoles fixed to their fiducial values. For example, in the t = 4 panel, the black curve is 
f(x) oc L(d\ . . . ,C^ E ,Cf E = x,C EE ,...). For comparison, naive pseudo-Q estimates are also 
shown with Gaussian errors (red curves). The pseudo-CV errors are noise only, while the condi- 
tional distributions include cosmic variance. 
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Fig. 8. — Conditional likelihoods for the I = 2 — 7 BB multipole moments (black curves), computed 
using the WMAP likelihood code by varying the multipole in question, with all other multipoles 
fixed to their fiducial values. For comparison, naive pseudo-Q estimates are also shown with Gaus- 
sian errors (red curves). The pseudo-C^ errors are noise only, while the conditional distributions 
include cosmic variance. Note the large difference between the likelihood code and the pseduo-C^ 
value for £ = 3; this mode is sensitive to the time-orderd data baseline and is extremely poorly 
measured by WMAP, illustrating the complicated noise structure of the polarization data on large 
scales. 
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Fig. 9. — WMAP 5-year EE power spectrum at low-£. The error bars are the 68% CL of the 
conditional likelihood of each multipole, with the other multipoles fixed at their fiducial theory 
values; the diamonds mark the peak of the conditional likelihood distribution. The error bars 
include noise and cosmic variance; t he point at I = 7 is t he 95% CL upper limit. The pink curve 
is the fiducial best-fit ACDM model (jPunklev et alJbood ). 



-26- 



0.006 
0.003 
i 0.002 [ 
— 0.001 

LU 
LU 



o 



o r - 

■0.001 - 
■0.002 



WMAP 5yr O 
Boomerang '03 O 
CAPMAP '08 
CBI '05 O 
DASI '05 O 
QUAD '07 O 




10 



100 500 

Multipole moment / 



1000 1500 



Fig. 10. — WMAP 5-year EE power s pectrum, compared with results fr om the Boomerang 



(jMontrov et alJ 12009. gre en). CAPMAP (IBischoff et al.1 120081. o range). CBI (jSievers et al.l 120071 . 
red), DASI (|Leitch et alJbood . blue), and QUAD (|p. Ade et al-lbooi purple) exp eriments . The 



pink curve is the best-fit theory spectrum from the ACDM/ WMAP Markov chain (jDunklev et al 



20081 ). Note that the y-axis is Cf h , not £(£ + l)Cf A /(2vr) 
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Fig. 11. — Left: One-dimensional marginalized likelihood distributions of erg for various treatments 
of the source uncertainty in the likelihood code: the standard likelihood function [black], the al- 
ternative treatment of the source uncertainty described in equation (|A5j) [blue], the alternative 
treatment, but with the unresolved point source error increased by x5 [cyan]. The agreement be- 
tween black and blue curves shows that the standard treatment is producing the correct answer. 
Right: ID marginalized likelihood distributions of n s for various treatments of the beam uncertain- 
ties: the standard likelihood function [black], the alternative treatment of the beam uncertainty 
described in equation (|A5P [red], the alternative treatment, but with the beam error increased by 
a factor of 20 [orange]. The agreement between the black and red curves shows that the standard 
treatment is producing the correct answer. 
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Fig. 12. — The unbinned WMAP 5-year temperature (TT) power spectrum (black), compared 
with the WMAP 3-year result (red). The slight upward shift of the 5-year spectrum relative to the 
3-year spectrum is due to the change in the beam transfer function. The pink curve is the best-fit 
ACDM model to the WMAP5 data. 
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Fig. 13. — The unbinned WMAP 5-year temperature (TT) power spectrum as a function of fre- 
quency, divided by the best-fit ACDM model to the WMAP data. 



